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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM NO. 1349 


ON A CLASS OF EXACT SOLUTIONS OF THE EQUATIONS OF MOTION 

OF A VISCOUS FLUID* 

By V. I. Yatseyev 


The general solution is obtained herein of the equations of motion 
of a viscous fluid in which the velocity field is inversely proportional 
to the distance from a certain point. Some particular cases of such 
motion are investigated. 

1. The motion of a viscous fluid with velocity field and pressure 
in spherical coordinates can be given by the following expressions: 
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A particular solution of the equations of Navier-Stokes for this case 
was obtained by Landau (reference l) . In the present paper a general 
solution is given of the equations of Navier-Stokes for the motion of 
the class under consideration. 


Substituting expressions (l) in the equations of Navier-Stokes and in 
the equation of continuity yields the following system: 


F 2 + f 2 - fF' + 2g 


ff' + g' - V 


+ F' cot 0 - 2f' - 2F - 2f cot 0 


V j^F" 

f" + f' cot 0 + 2F' - f (1 + cot 2 0) 


] = ° 


= 0 


F + f' + f cot 0 = 0 


( 2 ) 

(3) 

(4) 


Determining F from equation (4) and substituting in equations (2) and 
(3) give 
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f ,2 + ff .» + 3 ff . cot 0 + 2g - V + 2f" cot 0 - f ' (2 + cot 2 0) + 

f cot 0 (1 + cot 2 0)J = 0 (5) 

ff' + g' +V f" + f cot 0 - f (1 + cot 2 0)1 = 0 (6) 

Differentiating expression (6) 

f ' 2 + ff" + g" +V j^f" + f" cot 0 - 2f ' (l + cot 2 0) + 

2f cot 0 (1 + cot 2 0)J = 0 (7) 

Eliminating the nonlinear terms f' 2 , ff", and ff' from equation (5) 
with the aid of equations (6) and (7) yields a linear equation in the 
function g + 2Vf ' : 


(g + 2Vf * ) " + 3 cot 0 (g + 2Vf * ) ' - 2 (g + 2Vf ' ) = 0 


( 8 ) 


the general solution of which is in the form 


g + 2Vf> = 2P 2 l3 - cos - 9 (9) 
sin 0 

where 2P 2 a and 2P 2 b are constants of integration. 

Integrating equation (6) 

f 2 + 2g + 2V (f* + f cot 0) = - 2P 2 c (10) 

where 2V 2 c is the constant of integration. 

The function g($) is eliminated from equations (9) and (10) to 
give an equation of the Riccati type for the function f:^ 

( 11 ) 


= ^ f 2 + f cot 0 + 2V 




h cos 0 - a c 
sin 2 0 2 


^After sending the manuscript to press the author obtained from 
L. D. Landau a communication on the work of N. Slezkin (reference 2) in 
which he arrived at the same equation by a different method. 
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The substitution 


f = - 2VX'(0)/X(e) 
reduces equation (ll) to the linear equation: 


( 12 ) 


X” - X' 


cot 9 + 


( b cos 9 - a 
\ sin 2 9 


+ 



(13) 


which by the substitution 


z = cos^ (0/2) 

is transformed into an equation of the Fuchs ian type: 


(14) 


d 2 X a + b - 2 (b + c) z + 2cz 2 
dz 2 4z 2 (z - l) 


(15) 


The usual computations (reference 3), which are omitted herein, give the 
general solution of equation (15) as: 


x(e) 


( 9\r ( . e\i+o*3-r f / 

(cos g] Sinj <cjF I a,|3,r> 

C2F ( a + 1 - y, & + 1 - y,2 - x> cos 2 


cos*- 




(16) 


where the parameters of the hypergeometric function a, (which can 
also have complex values) are connected with the constants of integra- 
tion a,b,c by the formulas: 


a =T 2 - (l + a + P) y + ^ i 

b = (a + 3 - 1) r - — ■ ■ +— - + | ► 

c = (a - P) 2 - 1 

2 J 


( 17 ) 
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Formulas (4), (9), (16), and (17) give the general solution, depend- 
ing on the four constants a, b, c, and A = C 2 / c^, of the Navier-Stokes 

equations for the class of motion of a viscous fluid under consideration. 
The constants of integration a, b, and c are expressed in terms of the 
corresponding tensor components of the density of the momentum transfer: 


^ v i ^ v k\ 

Hik = P8 ik + PVj.v k - pv [ — + —J 


Carrying out the computations 


n 


2v^p / b cos 9 - a 


W~ r 2 


sin^ 9 


_ _ 2v^p /a - b cos 9 c 

J ~ L (90 — •? I o o 

sin^ 9 d 
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The streamlines are determined by the equation: 


(18) 


(19) 


dr/v r = rde/v 0 (20) 

the integration of which gives 

const/r = f sin 9 (21) 

2. Attention is now given to two particular examples for which the 
equation of Fuchs degenerates. 

(a) Equation (15) has only one regular singular point, z = t*». In 
this case . 


a = b = c = 0 

and therefore by equations (19) 


( 22 ) 


(23) 


Hcpcp- - 0 
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The particular solution of equation (15) 


X(0) = 2z - 1 - A 


(24) 


leads by formulas ( 4 ), 


(9)^ and ( 12 ) 


to the solution found by Landau: 


F(0) = 2V 


A 2 - 1 

(A - cos 0 ) 2 


1 


f(e) 


2V sin 0 
cos 6 - A 


g( 0 ) = 4 V 2 1 _ - _ A cos 9 ( 25 ) 

(cos 0 - A ) 2 

This solution is analogous to the problem of a stream flowing out of the 
end of a thin pipe into a region filled with the same fluid. It is the 
only regular solution for all values of the angle 0 . 

(b) Equation (15) has only two regular points z = 0 and z =». 

In this case it follows from equation (15) that 


a = b = c / 0 (26) 

and equation (ll) becomes Euler's equation 


2 z 2 (d^/dz 2 ) - aX = 0 (27) 

the general solutions of which are 


Correspondingly, 

f(0): 


X(0) = e x / 2 cosh (nx + A) for a> - l/Z 
X( 0 ) = e x / 2 cos (nx + A) for a < - l/2 
X(0) = e x / 2 (1 + Ax) for a = - l /2 
the following equations are obtained for 


'N 

> (28) 
j 

the function . 


f = 2V rr 'cos 0 ( n tanh (^ + a) + 1/2J 

f = 2V TTTS I 0 fi ■ n tan (nx + A) ) 


f = 2V 


sin 0 j A n /A 
1 + cos 0 \l + Ax + X / 2 f 


for a > - l /2 
for a < - l /2 
for a = - l /2 


(29) 

(30) 

(31) 



6 


MCA TM 1349 


where x = In (l + cos 9 ) , n = i Jyi + 2a J. (For a = 0, n = l/2 in 
equation (29) the solution of Landau is again obtained.) 

For the solution of equation (29) by formula ( 4 ) 


F(0) = - 2V < n tanh (nx + A) + l/2 > + 2V 


1 - cos 9 
1 + cos 9 


n 


cosh (nx + A) 


(32) 


while g ( 9 ) is determined by formula (9) . 


The equation of the streamlines is in the form 


const/r = (l - cos 9 ) 

where the values of the constants 
conditions 


r 1 

<n tanh (nx + A) + l/2^ 
n and A are determined from 


(33). 

the 


f (§) = 2V (n tanh A + l/2) 




(34) 


The obtained solution corresponds to the problem of the stream flowing 
from the half line 9 = jt into a region filled with the same fluid. 

For solution (30), the parametric equation for the. . streamlines is 
in the form 


const/r = (2 - e x ) ^ !_/ 2 - n tan (nx + A) J 

0 = arc cos (e x - l) (35) 

The function (35) for 9 tt is a strongly oscillating one. It can 
therefore be concluded that solution (30) has no physical sense. 

The author acknowledges the suggestions and help received from 
Professor Y. B. Rumer. 


Translated by S. Reiss 

National Advisory Committee for Aeronautics 
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